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We describe a framework for designing novel materials, combining modern first-principles elec¬ 
tronic structure tools, materials databases, and evolutionary algorithms capable of exploring large 
configurational spaces. Guided by the chemical principles introduced by Antipov, et. al, for the 
design and synthesis of the Hg-based high-temperature superconductors, we apply our framework 
to design a new layered copper oxysulfide, Hg(CaS)2Cu02. We evaluate the prospects of supercon¬ 
ductivity in this oxysulfide using theories based on charge-transfer energies, orbital distillation and 
uniaxial strain. 


The superconductors with the highest known tran¬ 
sition temperatures at ambient pressure are all lay¬ 
ered compounds containing planes of copper and oxygen 
where the superconducting electrons reside, separated by 
“spacer layers” composed of other elements. For a given 
compound, varying the doping level to an optimal value 
near 0.15 holes per copper maximizes the superconduct¬ 
ing transition temperature T c . Since the copper oxide 
planes are a common ingredient, the large variabilty in 
optimal Tc s between compounds must then be controlled 
by the spacer layers, which function to tune the chem¬ 
ical and structural properties of the copper oxide layer. 
Finding novel compositions for the spacer layers is key to 
discovering new superconductors with higher transition 
temperatures. 

Theoretical design of new compounds is challenging 
due to the vast combinatorial space of elements and 
the large number of constraints: preferred oxidation 
state, electronegativity, atomic radii, preferred local co¬ 
ordination environment, and overall electrical neutrality. 
These microscopic properties in turn determine the local 
structural stability, global configurational minimum, and 
thermodynamic stability. Finding the global low-energy 
structure is the computational bottleneck. If a given 
composition results in a stable compound, we still need to 
determine whether it tunes the low-energy Hamiltonian 
so that T c is enhanced, which imposes further screening 
criteria. 

Rather than design a compound from scratch, we 
adopt a more moderate approach: we take as a starting 
point the family of cuprates with the highest transition 
temperatures, the Hg-based cuprates, and modulate its 
spacer layers. We benefit from the vast body of chemi¬ 
cal intuition accumulated for the cuprates, for instance, 
laid out especially clearly in Ref. 1, and use modern elec¬ 
tronic structure methods [2, 3], materials databases [4, 5] 
and evolutionary algorithms [6] to efficiently screen for 
compositions which have the desired structure and the 
best prospects for stability. We use three proposed theo¬ 
ries, based on epitaxial strain [7], the charge-transfer en¬ 
ergy [8] , and orbital distillation [9], to evaluate the resul¬ 
tant compound’s prospects for superconducivity. In this 
manuscript, we describe the general framework for guided 


design of new materials, and demonstrate its workflow 
by applying the principles to design a new layered cop¬ 
per oxysulfide Hg(CaS)2Cu02, which we abbreviate (HC- 
SCO). If it can be synthesized and doped, we believe 
HCSCO will be a high-temperature superconductor. 

RATIONAL DESIGN 

The cuprates are a functional stack: the composition 
of each layer is chosen to play a specific role (Fig. 1). 
The central copper oxide (CUO 2 ) plane supports super¬ 
conductivity and roughly constrains the in-plane lattice 
constant. The remaining layers must tune the chemical 
potential of the C 11 O 2 layer without rumpling the plane 
or introducing disorder, and isolate each CuC >2 plane to 
create a 2D system. 

In the Hg-cuprates, the HgOa layer harbors dopant 
atoms which tune the chemical potential. The BaO layers 
immediately adjacent to the C 11 O 2 plane spatially sepa¬ 
rate the superconducting electrons from the detrimental 
effects of the disordered dopant layer [10, 11]. Addition¬ 
ally, the highly ionic nature of the BaO layer means they 
do not capture dopant electrons intended for the C 11 O 2 
plane. The preference of Hg to be dumbbell coordinated 
bonds the entire structure together without introducing 
structural distortions [1]. Finally, the highly ionic O- 
Hg-0 dumbbells minimize c-axis hopping to maintain 2- 
dimensionality. 

Due to their spatial proximity, the adjacent BaO lay¬ 
ers tune the hoppings and interaction strengths of the 
in-plane Hamiltonian. Designing compounds with novel 
adjacent layers provides a mechanism for controlling su¬ 
perconductivity by, e.g., reducing the charge-transfer en¬ 
ergy. We quickly realized the most stringent constraint 
is structural stability, so we focused first on isolating 
stable candidates, then subsequently investigating their 
electronic properties. 

To maximize the likelihood that a proposed composi¬ 
tion is stable, we note that the layers adjacent to the 
Cu0 2 plane form a rock salt structure. Using materials 
databases, we selected all naturally occurring rock salt 
compounds AX, composed of an cation A and an an- 
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FIG. 1. Each layer in the structural stack of HgBa 2 Cu 04 -|_i performs a specific chemical, structural and electronic function. 
We focused on tuning the layers immediately adjacent to the CuCh plane which function to tune the in-plane hamiltonian, but 
must also be charge neutral and not interfere with doping. 


5.0 


Ba 


I_ 

1 



CiHg O 

_U__I_ 

2 3 4 

Electronegativity (Pauling) 


FIG. 2. Rock salt compositions AX plotted as a function of 
their lattice constant vs. electronegativity of the anion X, 
and color-coded by anion’s nominal valence. Electronegativi¬ 
ties of elements existing in HgBa 2 Cu 04 are marked. We seek 
compositions which have (a) lattice constants comparable to 
the in-plane Cu-Cu distance in the cuprates (~ 3.80-3.95 A), 
and (b) anion electronegativities greater than that of Cu, as 
dopants intended for the copper oxide plane would otherwise 
be captured by the anion, resulting in additional Fermi sur¬ 
faces. The compositions of interest are marked by a shaded 
rectangle. Arrows point to BaO and CaS. 


ion X (Fig. 2). The phase space is large and the rate 
limiting step is structural prediction, so we quickly pre¬ 
screen candidates by discarding compositions with (1) 
large lattice mismatches relative to the in-plane Cu-Cu 
distance, which we took to be 3.82 A, and (2) anions less 
electronegative than Cu, as these anions would capture 
dopants intended for the superconducting plane, produc¬ 
ing additional Fermi surfaces. 

Given the remaining list of potential compositions, we 
screened for stable compounds as follows: 

1. Local stability - substitute AX for BaO layer of 
HBCO structure forming Hg(AX)2Cu02, fully re¬ 
lax the structure, then compute the phonons. Dis¬ 


card if there are unstable phonons. 

2. Configurational minimum - place one formula unit 
of the elements (lHg, 2A, 2X, lCu 20) in a box 
and apply evolutionary algorithms to find most en¬ 
ergetically stable crystal structure. 

3. Thermodynamic stability - compute convex hull 
of Hg-A-X-Cu-0 system and determine whether 
Hg(AX) 2 Cu02 lies below the hull. Synthesis gen¬ 
erally proceeds in oxygen environments so the hull 
is computed with fixed Hg-A-X-Cu stoichiometry 
and varying oxygen concentration, parameterized 
by the chemical potential /i(O 2 ). 

Compositions which pass all these hurdles have a good 
chance of being synthesizble in bulk. However, these 
criteria aren’t rigid, as compositions which exhibit only 
weak phonon instabilities are likely synthesizable and 
possibly superconducting. Additionally, compositions 
that lie above the convex hull may also be synthesized 
via epitaxy or under high pressures. 


METHODS 

The pymatgen implementation of the Materials 
API [12] was used to access the Materials Project [4] 
database to select all structures with lattice constants 
within ~ 5% of 3.9A, the nominal Cu-Cu distance. 

Structural relaxations within DFT were performed 
with VASP [2, 3] using PAW potentials [13] and the 
PBEsol functional [14]. We computed phonons at the 
(0, 0), ( 7 r, 0) and (ir, 7 r) fc-points using 1x2 and \J2 x \/2 
supercells constructed with the help of ISOTROPY [15] 
to check local stability. 

We performed a fixed-composition evolutionary search 
for crystal structure using USPEX [16, 17], with VASP 
as the underlying DFT engine. A unit cell consisting of 
one formula unit and an initial population of 256 ran¬ 
dom structures were used, decreasing to 64 in following 
generations. USPEX successfully found existing cuprate 
structures (LSCO, YBCO, HBCO) in benchmark tests, 
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demonstrating its efficacy in the few-atom (8 atoms/f.u.) 
but many-species (5 elements) regime. 

The Gibbs phase diagrams for the Hg-A-X-Cu system 
at varying oxygen potentials / 4 O 2 ) were constructed us¬ 
ing the phase diagram implementation in pymatgen [18]. 
We cross-checked these results with the PBE functional, 
which is known to give improved atomization ener¬ 
gies [14], finding no significant differences in the relative 
energies between compounds. 

Bandstructure orbital character was analyzed with 
WIEN2k [19] and structures plotted using VESTA [20]. 
Electronic parameters were extracted by downfolding to 
atomic-like orbitals, defined in [21]. 


RESULTS 


Our materials design framework efficiently isolates 
compositions with desired crystal structures and prop¬ 
erties. From roughly 350 rock salt compounds, pre¬ 
screening based on lattice constants and electronegativ¬ 
ity selected ~ 20 compositions. Screening for local sta¬ 
bility via phonon modes left three compositions: CaS, 
ZrAs, and YbS. Configurational stability using evolu¬ 
tionary structure prediction showed only a single com¬ 
position, Hg(CaS)2Cu02, adopted the HBCO structure. 
This compound has sulfur instead of oxygen as the api¬ 
cal atom and has the smaller calcium ion instead barium 
in the adjacent layer. We also performed structure pre¬ 
diction with two formula units, finding the same HBCO 
structure. Thermodynamic calculations show HCSCO 
is marginally unstable, lying 240 meV/atom above the 
convex hull at /r(02) = —16.48 eV, with an expected de¬ 
composition into Hg, Cu, CaS, and oxygen gas (Fig. 3). 

Turning to electronic structure (Fig. 4), we observe a 
single band crossing the Fermi level with strong d x 2_ y 2 
character, similar to other single-layer cuprates. The 
Fermi surface extremely 2-dimensional, and comparable 
to HBCO. Table I summarizes a comparison of the salient 
electronic parameters extracted using downfolding be¬ 
tween HCSCO and the single-layer HBCO superconduc¬ 
tor. 


Structure (A) 

U C d a pi ca l 

HCSCO 3.833 10.054 2.746 
HBCO 3.879 9.516 2.789 


Parameters (eV) 

A pc t A E z tpd tpp tpp 

1.94 0.20 1.35 0.61 0.09 

1.93 1.87 1.25 0.65 0.16 


TABLE I. Comparison of structural and electronic parameters 
of HgBa 2 Cu0 4 and Hg(CaS) 2 Cu0 2 . 



FIG. 3. The Gibbs phase diagram of the Hg-Ca-S-Cu system 
at an oxygen chemical potential of /r(0 2 ) = —16.48 eV, chosen 
because Hg(CaS) 2 Cu0 2 (red dot) is least unstable at this 
value. The phase diagram forms a tetrahedron with S, Hg, 
Cu and CaO at the vertices (elemental Ca isn’t stable under 
this oxygen environment). HCSCO lies in the interior of the 
tetrahedron, on the triangular face formed by Hg, Cu and 
CaS. 


DISCUSSION 


We evaluate the prospects of superconductivity in HC¬ 
SCO based on three chemically-based proposals. Ap¬ 
plying in-plane compression induces large increases in 
transition temperatures, with dT c /d(\oga) ~ —600 K in 
LSCO [7] and —4000 K in HBCO [1], where a is the 
in-plane lattice constant. We are aware that many vari¬ 
ables are involved, especially that compression may in¬ 
duce buckling, which decreases T c [22]. However, assum¬ 
ing all other factors are held constant, the 1 . 6 % compres¬ 
sion of HCSCO relative to HBCO would cause a 60 K 
increase in T c . 

An alternative proposal argues that tuning the charge- 
transfer energy A p( j, which roughly corresponds to the 
effective Coulomb interaction, controls superconductiv¬ 
ity [ 8 ]. As shown in Table I, the negligible difference in 
A p d between HCSCO and HBCO, combined with simi¬ 
lar hopping amplitudes, implies comparable T c ’s (Fig. 5). 
We expected that the increased polarizability of sulfur as 
compared to oxygen would reduce A p d, an effect which 
may be captured in methods beyond LDA, like the GW 
approximation. 

Finally, the orbital distillation proposal argues that a 
large admixture of the apical orbitals, in particular the 
Cu-d z 2 orbital, into the d x 2_ y 2 band suppresses T c [9]. 
Although the energy splitting A E z = E x <s _ y 2 — E z 2 is 
quite small (Table I), there is no admixture of apical or¬ 
bital charcter into the in-plane band (Fig. 4). Orbital 
distillation predicts T c ’s comparable to the HBCO su¬ 
perconductors. 

Turning to stability against phase separation, our cal- 
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FIG. 4. Bandstructure of HCSCO, with d x 2_ y 2 (top) and d z 2 
(bottom) character highlighted. A single band crosses the 
Fermi level of mainly d x 2_ y 2 character with little admixture 
of d z 2 character. 


culations indicate HCSCO lies slightly above the convex 
hull, well within the systematic errors of the theory used 
to estimate total energies. Furthermore it is known that 
many interesting functional materials are metastable [25] , 
protected from decay by large energetic barriers. In solid 
state synthesis, the Hg-cuprates are formed in the pres¬ 
ence of gaseous Hg [26] in addition to 0 2 , providing an 
additional tuning parameter which may enhance stabil¬ 
ity. Finally, prior experience with materials design indi¬ 
cates that compounds which are too stable may not be 
easily dopable [27, 28]. These considerations suggest av¬ 
enues for synthesis of HCSCO should be explored, both 
under pressure in bulk and especially via molecular beam 
epitaxy. 



(eV) 

FIG. 5. Scatter-plot of the maximal superconducting tran¬ 
sition temperature T C}lnax vs. charge-transfer energy A p d for 
the copper oxides (see Ref. 8 for details). The vertical bar 
marks the of HCSCO, which is close to that of the other 
mercury cuprates (labeled 1, 2 and 3 by the number of Cu0 2 
layers); thus, we expect similar transition temperatures. For 
comparison, we have also plotted the charge-transfer energies 
of two layered nickelates, LasNioOg [23] and LaiNHOs [24], 
where the Ni valence can be tuned near 3d 9 and are predicted 
to be nonsuperconducting (red circles). 


CONCLUSIONS 


We have created a workflow for systematically de¬ 
signing novel compositions based on existing compound 
families. Guided by chemical intuition and materials 
databases, our workflow quickly isolates promising com¬ 
positions for in-depth examination. Combining electronic 
structure tools, evolutionary algorithms, and methods for 
constructing phase diagram, we screen for local, configu¬ 
rational and thermodynamic stability to isolate new com¬ 
pounds with a high probability of synthesis. 

We have applied this workflow to design a new layered 
copper oxysulfide, Hg(CaS) 2 Cu0 2 , which we believe will 
be a high-temperature superconductor, with T c ’s compa¬ 
rable to those of the mercury cuprates. In the future, we 
expect the multi-layer family members of HCSCO will be 
designed, with formulas Hg(CaS) 2 Rn-i(Cu0 2 ) n where R 
is a 2+ cation placed between the Cu0 2 planes in the 
“infinite-layer” stack, as well as systematic extensions to 
the compositions of other cuprate families. 

To date, cuprate superconductors exist with oxygen, 
chlorine and fluorine in the apical position. The synthesis 
HCSCO would produce the exciting prospect of adding 
a fourth, sulfur-based family. 
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